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Abstract. We review our attempts to discover lost baryons at low redshift with "X-ray forest" of 
absorption lines from the warm-hot intergalactic medium. We discuss the best evidence to date 
along the Mrk 421 sightline. We then discuss the missing baryons in the Local Group and the 
significance of the z=0 absorption systems in X-ray spectra. We argue that the debate over the 
Galactic vs. extragalactic origin of the z=0 systems is premature as these systems likely contain 
both components. Observations with next generation X-ray missions such as Constellation-X 
and XEUS will be crucial to map out the warm-hot intergalactic medium. 
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1. Introduction: The Missing Baryons 

Given the title of this symposium, it should not come as a surprise to know that 
most of the baryons at low redshift are missing; here is how we know. The theory of 
big-bang nucleosynthesis, together with the observations of the deuterium abundance 
implies that the ratio of baryon density to the critical density of the Universe is about 
4%: 57b ~ 0.04/iyQ^ (Buries & Tytler 1998). The concordance cosmology with the new 
WMAP observations also predicts a similar number for fib (Bennett et al. 2003). This 
theoretical expectation is consistent with the observations of Lyman a forest at high 
redshift. At low redshift, however, the Lyman a forest thins out, so the baryons clearly 
do not stay in the warm, photoionized intergalactic medium at low redshift, as they did 
at high redshift. Perhaps, as the Universe evolved, and more and more galaxies formed, 
the baryons were accreted onto the virialized systems. The baryon census at low redshift, 
however, proved it not to be the case. All the stars and gas in galaxies account for only 
10% of all the baryons (Fukugita, Hogan & Peebles 1998). Hot gas in clusters of galaxies 
makes a similar contribution to Vlb- Cold gas at low redshift also adds to the baryon 
budget, but most of the baryons were found to be lost. This is the so called missing 
baryon problem. 

Cosmological hydrodynamic simulations found that as the Universe evolves, the gas in 
the intergalactic medium gets shock heated to high temperatures (10^~10^K) and forms 
a web like large-scale structure of filaments and sheets. Galaxies and clusters of galaxies 
form in dense knots in the filamentary structure, but even at low redshift most of the 
baryons still reside in the diffuse intergalactic medium (IGM). The temperature of the 
IGM is too high to keep hydrogen neutral, which is why Lya forest thins out. Given the 
predicted temperature range, the low redshift IGM was called warm-hot intergalactic 
medium, or WHIM, by Gen & Ostriker (1999), which is hotter than the warm IGM at 
high redshift, but cooler than the hot gas in clusters of galaxies. 

Even though most of the hydrogen in the WHIM is ionized, most of the metals are not 
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Figure 1. The fraction of oxygen in different ionization states vs. temperature. The top left 
panel is for the pure collisional ionization while photoionization from metagalactic UV and X-ray 
background is included in the other panels in which the density is lO""*, 10~^, and 10~®cm~^ 
from top right to bottom right. 

fully ionized. Thus, in spectra of distant quasars, the WHIM may imprint its signature by 
absorption lines of highly ionized metals. Oxygen being one of the most abundant metals, 
absorption lines of Ovii arc expected to be strongest (see fig.l, adapted from Mathur 
et al. 2003. O VII is the dominant ionization state of oxygen in the WHIM temperature 
range). 

2. The Quest for WHIM 

The resonance transitions of the dominant ions fall in the X-ray band-pass. The 
strongest line of O vii , for example, is at 21.602A. One thus expects to detect an "X-ray 
forest" of absorption lines from hydrogen-like and helium-like ions in X-ray spectra of 
distant quasars. Such an experiment, however, could not be performed before the launch 
of Chandra and XMM-Newton observatories, which, for the first time, carry gratings on 
board allowing high resolution spectroscopy. Even then, this is a very difficult experiment 
because these X-ray telescopes are small and given how faint the quasars are, the S/N 
of spectra becomes a limiting factor in detecting weak and narrow absorption lines due 
to WHIM. For this reason we came up with a strategy of observing blazars in their high 
state, which would provide temporarily bright light beacons. 

2.1. Chandra, FUSE and HST observations of Mrk 421 

One siicli blazar is Mrk 421 which we monitored with the Rossi X-ray Timing Explorer 
(RXTE). When the Mrk 421 flux was almost two orders of magnitude larger than its 
quiescent flux, we triggered a target-of-opportunity (TOO) observation with Chandra low 
energy transmission grating (LETG). Thus we obtained the highest S/N grating spectrum 
among all targets ever observed with Chandra (see figure 2, adapted from Nicastro et al. 
2005a, b). The black dots with error bars are the data and the blue line running through 
the spectrum is the Mrk 421 continuum folded with the instrumental response. The 
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Figure 2. The Chandra LETG spectrum of Mrk 421. Note the detection of absorption lines 
from various H-like and He-line ions at redshift zero and at z=0.011 and z=0.027 marked with 
vertical bars. 



red lines mark the absorption lines in the spectrum which we identify with transitions 
of several ions at redshifts zero, 0.011 and 0.027. The ions are labeled and the three 
redshifts are marked with vertical bars. The middle left panel shows the strong transition 
of Ovil Ka which is clearly detected at all the three redshifts. The strongest line is at 
z=0, and we will come to this in a later section. 

We also observed Mrk 421 with FUSE to look for any possible absorption lines of O vi. 
The FUSE spectrum, shown in figure 3 (adapted from Nicastro et al. 2005b), is rich with 
absorption lines from the interstellar medium of the Galaxy. However, O vi absorption 
lines from the two intervening X-ray absorption systems are not detected. 

Mrk 421 was also observed with HST and a portion of the STIS spectrum is shown in 
figure 4 (adapted from Nicastro et al. 2005b). Again, we do not detect HI Lya absorption 
from the two intervening X-ray systems. 

To summarize, we detect two intervening absorption systems toward Mrk 421. The 
z=0.011 system is detected through strong absorption lines of Ovii and Nvii , though 
O VI and H i are not detected. The z=0.027 system is detected through strong absorption 
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Figure 3. The FUSE spectrum of Mrk 421. No O vi absorption lines are detected at the 
redshifts of intervening X-ray absorption systems. 



lines of Ovii, Cvi, Nvi, Nvii and less significant lines of Neix and Ovii K/3. Ovi is 
not detected in this systems as well and Hi detection is also not strong, if at all present. 

With the column density of detected ions and the upper limits obtained from non- 
detections we determine the physical properties of the absorbing gas. The method can be 
understood from figure 1 in which fractional ionization of oxygen is plotted as a function 
of temperature. The top left panel corresponds to pure coUisional ionization; we sec that 
the fraction of oxygen in a given ionization state is a strong function of temperature 
and higher ionization states reach their peak fraction at higher temperatures. Thus, the 
column density ratios of two ions depend upon temperature. 

In the other three panels of figure 1, photoionization from metagalactic UV and X-ray 
backgrounds is also included in calculating ionization fractions. As expected, the shapes 
of the curves change when photoionization is added. The three panels correspond to three 
different densities of the IGM, which are affected by photoionization to different degrees. 
Thus, the ionic ratios become not only dependent on the temperature, but also on the 
density of the IGM. In particular, we find that the O vii/0 vi ratio is a sensitive probe of 
the temperature of the IGM, while the Oviii/Ovi ratio measures the density (Mathur 
et al. 2003). This can be understood intuitively; if a system contains both Oviii and 
Ovi with measurable amounts, then photoionization must play a role, and lower the 
density, higher the effect of the photoionization. 

We can then make figures similar to figure 1 for other elements and generate expected 
values of column density ratios for a variety of ions. We can then compare our observed 
line ratios with the theoretical calculations to constrain the properties of the intervening 
systems. One such diagnostic plot is shown in figure 5 (adapted from Nicastro et al. 
2005b) in which various ion ratios are plotted vs. temperature. The thin curves represent 
the theoretical models while the thick curves represent observed values within la errors. 
All the observations overlap in the temperature range of 0.4-3.3 x lO^K. 

Strong detection of oxygen and non-detection of H I allow us to put a constraint on 
the abundance ratio to be [O/H] > — 1.46, a reasonable value for the WHIM. The total 
equivalent column density is Nh = 0.3-7.3 x 10^^ x lOl'-'/^l-i cm~^. We have no direct 
constraint on the density, but if it of the order of 10~^cm~^, then the inferred path length 
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Figure 4. The HST spectrum of Mrk 421 showing no HI absorption from the two intervening 

X-ray systems 

is of the order of a Mpc: D=0. 1-2.4 x 10['^/^l-i(n_5)~^. Thus all the inferred parameters 
of the absorption systems are consistent with those of the WHIM. 

We have detected two intervening O vii systems in the pathlength out to Mrk 421. This 
is just one sighthne and there are only two systems in it, so we are in the regime of low 
number statistics, but nonetheless we can calculate the number density of O vii filaments, 
albeit with a large error bar. In figure 6 we compare this with the theoretical curve and 
again find consistency. The theoretical curve in new simulations (Cen & Fang 2006) has 
lower normalization, but the Mrk 421 data point is still consistent with models within 
2a. 

2.2. Lost is Found? 

Given the consistency of all the observed parameters with the theoretical expectations 
for the WHIM, it is reasonable to assume that the observed systems are indeed WHIM 
filaments. We can then calculate the baryon density associated with these systems and 
compare it with the missing baryons. We find 

- o.o32l°:^^2 ^ ^Q-[o/H]., 

consistent with the missing flj,. We are, however, cautious about claiming that the lost 
baryons have been found because, once again, the errors on our number are large due 
to small number statistics. It is therefore imperative to observe more sightlines and 
detect more WHIM systems. We observed a sightline toward 1ES1028 + 511 with XMM- 
Newton and have tentative detections of a couple of more systems. 1ES1028 + 511 is at a 
higher redshift (z=0.361) than Mrk 421 (z=0.03), resulting in a smaller number density 
of O VII systems along this sightline. Note, however, that the column density threshold in 
the XMM-Newton spectrum of 1ES1028 + 511 is much larger and the observed number 
density is still consistent with the theoretical curve in figure 6. If these are confirmed, 
then our solitary point on the curve in figure 6 will have company and the errors on fif, 
will be lowered. 

3. The z=0 Absorption Systems 

In high quality spectra of several extragalactic sources, absorption lines of O vii Ka are 
detected at redshift zero (see figure 2 for the z=0 O vii line in the Mrk 421 sightline). The 
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Figure 5. Constraints on tlie z=0.011 system from observed line ratios. The thin lines are the 
theoretical models while the thick portions on the lines mark observed values with errors. The 
temperature range where all observations overlap is around lO^'K. 

origin of such absorption systems is debated in the literature; is it from the ISM or the halo 
of our own galaxy or is it from larger scale structures such as the Local Group? It should 
be noted that O vi absorption is also ubiquitous in extragalactic sightlines observed by 
FUSE. One normally detects a broad component of O VI at redshift zero from the thick 
disk of our own Galaxy and often an additional component at a velocity offset of more 
than 100 km/s from the so called high velocity cloud. The spectral resolution in the X-ray 
band is not high enough to resolve the components seen in the FUV. It is not immediately 
obvious, therefore, whether the O vii line is related to any of the O vi components. The 
answer to this question is related to the origin of O vii systems. 

3.1. Galactic or Local Group? 

This has become such an important question because we would like to know the associated 
mass. If the z—0 systems are from the local large scale structure, then the associated mass 
can be huge. It may even account for the missing baryonic mass in the Local Group. On 
the other hand, if the z~Q O vii systems are from the halo of our Galaxy, then they may 
carry insignificant amount of mass. 

Theoretically, some models of galaxy formation and evolution predict hot halos around 
galaxies. Similarly, as mentioned in §1, in the cosmological simulations galaxies form 
in dense knots in IGM filaments, so we expect to find a local WHIM filament(s) in 
which our Galaxy is embedded. This expectation was further confirmed with constrained 
simulations of the local volume (Kravtsov et al. 2002) which show large reservoirs of 
O VII gas in extended structures around the Galaxy. Thus, understanding the origin of 
O VII systems will help test the models of galaxy formation and/or large scale structure. 

We have been investigating this question systematically (Nicastro et al. 2002, Williams 
et al. 2005, 2G06b, 2007) by modehng the observed Ovii and Ovi systems. We find 
interesting differences among different sightlines implying complexity in the structure 
surrounding the Galaxy. We could not determine the origin of the O vii systems; neither 
Galactic nor extragalactic origins could be ruled out. One conclusion that we could reach 
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Figure 6. The cumulative number density of Ovil filaments due to WHIM is plotted as a 
function of observational column density threshold. The theoretical curve is from Fang et al. 
2002 and the data point with error bars is from our observations. 

is that the Ovii absorption Unes do not arise in the gas producing the Ovi systems. 
This is most directly seen in the Mrk 279 sightline (figure 7, adapted from Williams et 
al. 2006) where we found a velocity offset between the O Vii line and the Ovi HVC. 

To summarize, Ovil absorption is not associated with any single Ovi component. 
Neither Galactic or extragalactic origin of the O vii systems can be ruled out. Most 
likely, a variety of phenomena are responsible for the O Vii systems as indicated by the 
differences in their "b" values (velocity dispersion parameters). Better, higher resolution, 
data are crucial in understanding the origin of O VII systems. 



4. Conclusions & Future Prospects 

As mentioned by C. Frenk (these proceedings) the claims of the discovery of WHIM 
have been controversial. The best observed evidence, nonetheless, is from the Mrk 421 
sightline in which the chance probability of detecting the observed absorption features 
by statistical fluke is less than 1%. Detections with higher significance are desirable, of 
course. Given the importance of the WHIM discovery, we need to observe more sigh-lines 
and solidly detect at least a few more X-ray forest systems before counting the baryons. 
The community needs to support observing programs with these goals. It is worth noting 
at this point that because of the viewing angle constraints on Chandra observations, TOO 
observations of blazars in high state is not a viable strategy any more. 

Even if we succeed in our effort to detect WHIM with Chandra or XMM-Newton , 
we will be only observing the tip of the iceberg (see figure 8). To map out the WHIM 
filaments at low redshift, by going down to the column density thresholds much more 
representative of the large scale structure, Constellation-X and XEUS will be essential. 
These are future missions of NASA and ESA respectively; mapping the low redshift 
WHIM is one of their major science drivers. We need to obtain a definite proof of concept 
with Chandra and XMM-Newton before these mission concepts are finalized. 

The science of the z=0 Ovil systems is equally interesting which bears directly on 
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Figure 7. The X-ray and UV z=0 systems toward Mrk 279 show velocity offset 

the models of galaxy formation and/or models of large scale structures. The spectral 
resolution offered by Chandra and XMM-Newton gratings is just not good enough to 
imdcrstand the origin of these systems. The debate over the Galactic vs. cxtragalactic 
origin of O Vii systems is premature because we do not have data quality to resolve the 
issue. Given what we understand from other poor groups similar to the Local Group, the 
O VII systems must contain both of these components. We propose that the O vii lines 
contain two (or more) components (Mathur et al. 2007) which will be resolved in higher 
resolution spectra. Future higher resolution observations are crucial to test this prediction 
and to make significant progress in this field. 
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Figure 8. The O vil column density map from tiie cosmological iiydrodynamic simulation of the 
WHIM (from Chen et al. 2003). The top left panel has a column density threshold of 10^* cm~^. 
The other panels from top right to bottom right correspond roughly to detection thresholds of 
XEUS, Constellation-X, and Chandra / XMM-Newton respectively (assuming metallicity of a 
tenth solar). 

Discussion 

BOMANS: The intervening absorption lines in the Chandra spectrum of Mrk 421 are not 
detected with XMM-Newton . Could you comment on that? 

Mathur: Because of the higher S/N in XMM-Newton spectrum, it was natural to expect 
that the Chandra lines will be detected in the XMM-Newton spectrum. However, S/N 
is not the only factor in deciding the line detection threshold. The shape of the line 
response function makes a big difference, just the same way a point spread function 
affects the point source detection in imaging observations. The line spread function of 
XMM-Newton gratings has wings, increasing the equivalent width limit of line detections. 
There are a couple of more factors, which we have detailed in Williams et al. 2006a. The 
bottom line is that XMM-Newton non-detections are consistent with Chandra detections. 
Later, Rasmussen et al. re-analyzed the XMM-Newton data. Their upper limits are again 
consistent with Chandra detections, given the values in the paper, even though the text 
of their paper claims inconsistency! 

SCHAYE: We do not know the metallicity of the WHIM, so we can never know whether 
the missing baryons are found. It would be much more fruitful to assume that the missing 
baryons are found and use that to constrain the UV and X-ray background 

Mathur: Good point! However, the errors on the observations at present are too large 
to constrain the background parameters. Secondly, if we actually detect hydrogen with 
UV spectroscopy from the WHIM systems, then we can measure metallicity. 



